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The birth of stars involves not only accretion 
but also, counter-intuitively, the expulsion of mat- 
ter in the form of highly supersonic outflows^ -^. 
Although this phenomenon has been seen in young 
stars, a fundamental question is whether it also oc- 
curs amongst newborn brown dwarfs: these are the 
so-called 'failed stars', with masses between stars 
and planets, that never manage to reach temper- 
atures high enough for normal hydrogen fusion to 
occur'^. Recently, evidence for accretion in young 
brown dwarfs has mounted^^, and their spectra 
show lines that are suggestive of outflows^"^. Here 
we report spectro-astrometric data that spatially 
resolve an outflow from a brown dwarf. The out- 
flow's characteristics appear similar to, but on a 
smaller scale than, outflows from normal young 
stars. This result suggests that the outflow mech- 
anism is universal, and perhaps relevant even to 
the formation of planets. 

The nearby p-Ophiuchi Cloud (at a distance of 125 
pc; ref. 10) is an excellent example of a stellar nursery. It 
was as part of a near-infrared survey of this region^ ^ that 
p-Oph 102 was first detected. Its mass is given as 60Mj 
(ref.7), where IMj equals the mass of Jupiter or O.OOIMq. 
This places it firmly within the brown dwarf mass range 
(O.O13M0 < M < O.O75M0)3 and indeed this object 
has been spectroscopically confirmed to be a BD^. There 
is also strong evidence for the presence of an accretion 
disk^ and, in fact, the accretion rate (derived from its Ha 
profile)^ is estimated to be approximately 10^^ AIq yr^^. 
Finally it has been noted that its spectrum contains a 
number of forbidden emission lines suggesting an outflow^. 

Tracing an outflow from a young star is done using a 
variety of techniques depending upon wavelength. When 
the star itself is sufficiently evolved to be optically visi- 
ble, i.e. much of the surrounding natal gas and dust has 
been driven away, the outflow can be followed almost right 
back to the star through its permitted and forbidden line 
emission^. Nearest the source, the light from the coUi- 
mated outflow is most intense, giving rise to the so-called 
"micro-jets" observed, for example, from T Tauri stars, 
the precursors to stars like our sun^^. In extreme cases, the 
outflow is observed only very close (<1") to the source^^, 
i.e. within the typical "seeing" disc for ground based tele- 
scopes. It is under such circumstances that the technique 
of spectro-astrometry comes to the fore (see Methods and 
Supplementary Information). 



Our high resolution echelle observations (see Methods 
and Fig. 1) of p-Oph 102 first of all show its forbidden 
emission lines, e.g. [01] AA6300,6363 and [SII]A6731, are 
moderately blueshifted (V^ ~ -45kms~^). Such radial ve- 
locities are very similar to those seen in outflows from clas- 
sical T Tauri stars, and might be expected of an outflow 
from a young BD, as its escape velocity is similar^*. We 
note also the absence of a redshifted outflow component. 
This is again typical of T Tauri stars and is interpreted in 
terms of an obscuring disk that hides the redshifted com- 
ponent from our view^'^. We can thus immediately infer 
the presence of a disk around this brown dwarf (see below 
for a limit on its size). Such a conclusion is also in line 
with the mid-infrared excess seen from this object"*. Note 
that all velocities are with respect to the systemic veloc- 
ity of the BD. The latter was derived from the Li A6708 
photospheric absorption line and equals 7±8 kms^^ in the 
Local Standard of Rest frame. 

Another indication that an outflow is present comes 
from a cursory examination of the Ha line (Fig. 2). Its 
profile is clearly asymmetrical, i.e. the blueshifted wing of 
the line appears to be absorbed in a P-Cygni like fashion. 
We note however that we do not observe a classical P- 
Cygni profile, i.e one that dips below the continuum. Such 
a profile is in any event a rare occurrence even amongst T 
Tauri stars^^. 

If we are dealing with a scaled-down version of the 
outflow phenomenon seen in T Tauri stars then we expect 
the centre of emission in forbidden lines, to be spatially 
offset from the continuum, i.e. the BD. This offset is due 
to the fact that, for a coUimated outflow, such lines are 
quenched close to the star once the electron density be- 
comes high enough*^. In the case of outflows from T Tauri 
stars, typical offsets of 30-75 AU (0'.'2-0'.'5 at 150 parsecs) 
are seen, for example, in [OI]AA6300,6363 and the red [SII] 
doublet*^. If we assume brown dwarf outflows have similar 
opening angles and velocities to those from T Tauri stars, 
then the point at which the critical density is reached, 
naively scales with M°J where Mjet is the jet's mass loss 
rate (see Supplementary Information). Assuming the lat- 
ter depends linearly on the accretion rate, we would expect 
typical spatial offsets to be 3-10 times smaller in brown 
dwarf outflows in comparison to those from T Tauri stars. 

Spectro-astrometric (emission centroid offset versus 
velocity) plots are shown in Fig. 1 for the [OIJA6300, 
[OIJA6363, Ha and [SIIJA6731 lines. Because of high elec- 
tron densities close to the BD, [SIIJA6717 was too weak 
to provide a usable spectro-astrometric signal. Here the 
spatial offsets were measured after continuum subtraction 
(see the Methods and Supplemenatry Information for de- 
tails). 

We now consider the main results from various offset 
versus velocity plots. First, the centroids of all the mea- 
surable forbidden lines are displaced to the south, i.e. have 
negative offsets with respect to the continuum. These off- 
sets reach a maximum, of O'.'OS— O'.'l, at a blue-shifted ve- 
locity of approximately -40kms^^. We have already noted 
the absence of any corresponding redshifted emission and 
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that this suggests the presence of an obscuring disk. The 
scale of the blueshifted offset would suggest a minimum 
(projected) disk radius of O'.'l (> 15 AU at the distance 
of the p Ophiuchi Cloud) in order to hide any redshifted 
component. 

Second, there is no clear spatial offset in Ha even 
though its higher signal to noise potentially allows us to 
measure even smaller offsets than observed in the forbid- 
den lines. This is in agreement with the idea that most of 
the Ha emission arises from accretion^ on much smaller 
scales than are being probed here. It is also worth pointing 
out that if, as is almost certainly the case, the blue-ward 
dip in the Ha profile is caused by P Cygni-like absorp- 
tion from an outflowing wind in front of the star, no off'set 
should be expected. 

Third, both the line profiles and the spectro- 
astrometric signatures are very similar (albeit on some- 
what smaller scales) to what is seen in T Tauri stars with 
"micro-jets" . In particular the observed velocities and off- 
sets in the various forbidden lines are within the range 
we would expect for a coUimated outflow from a brown 
dwarf (see also Supplementary Information). We suggest 
that direct imaging, (using, for example, Fabry-Perot sys- 
tems) of this and other candidate brown dwarf outflows, 
should now be attempted. Such observations however will 
be very challenging, even with large telescopes, because of 
the expected faintness of the outflow. 

Methods 

Spectro-astrometry 

Conceptually the principles of spectro-astrometry are easy 
to understand. The profile of a star is smeared by atmo- 
spheric turbulence to appear gaussian (at least to a first 
approximation) rather than point-like. Whereas the width 
of the profile is determined by the so-called seeing, how 
accurately we can determine the centroicl of emission is, in 
theory for fixed seeing, limited only by the strength of the 
observed signal to noise ratio. Increasing the total number 
of detected photons increases the positional, or astromet- 
ric, accuracy, so that, in principle, milliarcsecond precision 
is possible with very large ground based telescopes^^"^^. 

Consider now a long-slit spectrum of a close binary 
system consisting of two virtually identical stars. We will 
assume that the slit is orientated along the same position 
angle as the binary. (We note that strictly this is not nec- 
essary: it is only necessary that the slit is not orthogonal). 
If the separation of the binary is considerably less than 
the seeing, the profile of the system in the spatial direc- 
tion will consist of a single gaussian-that is, the system is 
unresolved and the centroid of emission will lie exactly be- 
tween the two components. Now suppose that one of the 
two stars differs slightly from the other in being a strong 
Ha emitter; in such a case, the emission centroid will shift 
towards that star in the spectrum at the position of the 
Ha line. In this way it is possible to resolve certain types 
of binaries with separations well within the seeing limit^^ . 



In the case of a jet (pure emission line region) plus star 
(continuum source), one can go further and interpolate the 
continuum across a line, thereby allowing its contribution 
to be removed. It is then possible to measure separately 
the spatial centroid of the pure emission line region and 
determine its offset with respect to the continuum, that 
is, the parent star. Moreover, as the line can be emitted 
over a range of wavelengths, owing to the Doppler effect, 
it may also be feasible to recover spatio-kinematic infor- 
mation. For example, if the jet is bipolar, that is, it has 
oppositely directed blue- and redshifted flows from the 
source; the emission centroid of the red and blue wings of 
the line will be displaced to opposite sides of the contin- 
uum centre. 



The detailed method by which we measure offsets can 

briefly be described as follows. First, the centroid of the 
continuum emission in the spatial direction is determined 
using a one-dimensional gaussian fit. The line of such cen- 
troids, in the dispersion direction but excluding any re- 
gion where emission lines are present, is then fitted with 
a second-order polynomial, over a range of typically 200- 
300. In this way, instrumental curvature and tilting, with a 
characteristic frequency many times larger than the width 
of any line, is determined. The fit, to the centre of the con- 
tinuum, is then subtracted from the actual measured cen- 
troids, leaving residuals that are evenly scattered about 
the abscissa (that is, the fit defines the zero offset line). 
The continuum data points shown in Fig. 1 are thus the 
residuals. Finally, the two-dimensional fit to the contin- 
uum, broadened to take account of the point spread func- 
tion, is subtracted from the emission lines. Any emission 
line offsets are then measured. 



The accuracy (in arcseconcls) of the method is set by 
the error in the centroid of the gaussian fit, which depends 
on the seeing and the number of detected photons, A^. 
Formally, the error is given by Seeing/[2(2 In 2)f/2Nf/2], 
assuming that photon noise is the only source of noise. 
N, of course, is a function of the binning and the spatial 
sampling (pixel width). This explains why, for example, we 
can achieve a higher spectro-astrometric accuracy with a 
bright line, such as Ha than a weak one, for example, the 
[SIIJA6731 line. In some cases, it is necessary to bin up 
a weak line in the dispersion direction, as we have done 
to varying degrees for the [01] doublet and the [SIIJA6731 
line, to achieve sufficient signal to noise ratio. Note that 
we sometimes use different binning factors for the con- 
tinuum, in comparison with the line, so as to achieve a 
similar signal to noise ratio in both. This allows us to 
have comparable offset errors in both components, and 
to define the common la error lines shown in Fig. 1. As 
can be seen from the plots, the typical limiting offset that 
we can measure in the spatial direction (3tT) is around 
30mas. This corresponds to 4.5 au at the distance to the 
p Ophiuchi cloud. 
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Echelle spectroscopy 

The high resolution spectra of p Oph 102 were taken 
with the UV-visual Echelle Spectrograph (UVES) on the 
European Southern Observatory's 8m Kueycn Telescope, 
one of the telescopes in the Very Large Telescope (VLT) 
suite, in May 2003. A total of three 45niin exposures of 
the target were made, together with a series of flats and 
biases as well as an observation of an arc lamp for wave- 
length calibration. The slit was orientated north-south 
and had a width of 1 "while the seeing was 0'.'65. The 
central wavelength was set at 580nm, giving a spectral 
range of 450-680nm. Only the red part of the spectrum 
from 580 to 680nm, however, was analysed. The pixel 
scale was 0'.'182 and the spectral resolution R=40,000. The 
data were reduced using standard Image Reduction and 
Analysis Facility (IRAF) routines. 
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Fig. 1. Line profiles (top row) and spectro-astrometric plots (middle row) for the [OI]A6300 and [OIJA6363 doublet and 
spectro-astrometric plots (bottom row) for the Ha and [SIIJA6731 lines. The much narrower corresponding night-sky 
lines, peaking around -9kms~^, have been subtracted although no offset measurements were made in their vicinity 
(thus accounting for the data point gaps). Continuum and line offset points are represented by black triangles and 
squares respectively. All velocities are systemic and spatial offsets are in the north-south direction (in arcseconds) with 
negative offsets to the south. Dashed lines delineate the ilcr error envelope. For Ha , note the much smaller offset 
scale. The [SH] line is blue-shifted to around -40kms~"'^. 
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Fig. 2. The outflow signature in the Ha profile of p-Oph 102. The P Cygni-like dip in the hne profile is a strong 
signature of outflow activity. Ha emission from the brown dwarf is absorbed as it passes through material moving 
outwards along our line of sight. Because this material is moving towards us, the dip is on the blueward side of the 
line. Classical T Tauri stars are strong Ha emitters, and P Cygni Ha profiles originally confirmed that such protostars 
drive outflows. The dip in the Ha emission of p Oph 102 is at approximately the outflow radial velocity determined 
from the forbidden lines. 



